Draft version December 9, 2009 

Preprint typeset using 1^1^^ style cmulatcapj v. 03/07/07 



STRUCTURE AND MORPHOLOGIES OF Z - 7 - 8 GALAXIES FROM ULTRA-DEEP WFC3/IR IMAGING 

OF THE HUDF ^ 

P. A. Oesch^, R.J. BouwENS'\ CM. Carollo^, G.D. Illingworth^ , M. Trent:"', 
M. Stiavelli'"', D. Magee^, I. Labbe", M. Franx^ 

Draft version December 9, 2009 



o 
o 

(N 
o 

Q 



O 
U 

o 

cn 
> 
m 

00 

in 

0^ 
O 
0\ 
O 



X 



ABSTRACT 

We p resen t a first morpho l ogical study of z ~ 7 — 8 Lyman Break galaxies (LBGs) from lOesch et al.l 
(|2009bD and iBouwens et al.l (|2009af ) detected in ultra-deep near-infrared imaging of the Hubble Ultra 
Deep field (HUDF) by the HUDF09 program. With an average intrinsic size of 0.7 ± 0.3 kpc these 
galaxies are found to be extremely compact having an average observed surface brightness of /ij ~ 26 
mag arcsec"^, and only two out of the full sample of 16 z ~ 7 galaxies show extended features with 
resolved double cores. By comparison to lower redshift LBGs it is found that only little size evolution 
takes place from z ~ 7 to z ~ 6, while galaxies between z 4 — 5 show more extended wings in 
their apparent profiles. The average size scales as (1 -f z)^™ with m = 1.12 ± 0.17 for galaxies with 
luminosities in the range (0.3-l)L*^3 and with m = 1.32 ± 0.52 for (0.12-0.3)L*^3, consistent with 
galaxies having constant comoving sizes. The peak of the size distribution changes only slowly from 
z ~ 7 to z ~ 4. However, a tail of larger galaxies (> 1.2 kpc) is gradually built up towards later 
cosmic times, possibly via hierarchical build-up or via enhanced accretion of cold gas. Additionally, 
the average star-formation surface density of LBGs with luminosities (0.3-l)L*^3 is nearly constant at 
"^SFR = 1-9 M0yr~^kpc~^ over the entire redshift range z ~ 4 — 7 suggesting similar star-formation 
efficiencies at these early epochs. The above evolutionary trends seem to hold out to z ~ 8 though 
the sample is still small and possibly incomplete. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: structure 



1. INTRODUCTION 

The newly installed WFC3/IR camera on the Hub- 
ble Space Telescope has opened up a new territory in 
the study of galaxies at z > 6.5. Its increased ca- 
pability has lead to numerous detections of galaxies at 
z > 6.5 alr eady in the first epoch data of the HUDF09 
(lOesch et a l. 2009b; Bo uwens et al.ll2009aHMcLure et al.l 
l2009HBunker et al...2009h . 

Understanding the evolution of galaxy sizes and mor- 
phologies out to z > 6.5 can provide essential clues to 
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galaxy formation models. While hy drodynamical galax;v 
formation simulati ons of LBGs (e.g.lFinlator et al.lf2006l : 
iNight et all 120061: iNagamine et al.l l2008D have focussed 
on the prediction of the evolution of the galaxy luminos- 
ity and mass functions, the sizes of these galaxies have 
not been addressed in great detail so far. By providing 
these results we hope to stimulate interest in this key 
aspect of galaxies at early times. 

Based on semi-analytical modeling, the sizes of Lyman 
Break Galaxies (LBGs) at fixed luminosity are expected 
to slo wly decrease towards earlier cosmic times (Mo et al. 
[1998[ ). which is in good agreement with earlier observa- 
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tions of LBGs (e.g. iFerguson et al.|[200l : IBouwens et al.l 
[200l . 

The UV morphologies of LBGs have been studied be- 
tween z 2 — 6 from both ground-based and Hub- 
ble Space Telescope imaging. They have been found to 
be very compact, but often containin gj multiple cores, 
espec i ally at brighter magni t udes (e.g . Giavalisco et alj 

IRavindranath et al. 200^; 
20071 : IConselice fc ArnoTdl 
mechanism for star-formation 
in these galaxies is still debated. LBGs with multi- 
ple cores have been interpreted as merging systems, 
with star-formation triggered by this interaction, which 
is a natural scenario in a hierarchical un iverse (e.g. 
ISomerville et all 120011 : lOverzier et all l2008l ) . However, 
these individual cores could also originate from individ- 
ual star-forming clumps within a larger gas-dominated 
disk galaxy, whereby this clumpy state is kept alive due 
to g as streams replenis hing the disk galaxy with cold gas 
(e.g. iDekel et al.lf2009t) . 

Here we present the highest resolution observations 
available to date of galaxies from a time when the uni- 
verse was only ^ 800 Myr old. Our analysis is primarily 
based o n the sample of 16 z ~ 7 LBGs from lOesch et al.l 
()2009bf ) . We also discuss the 5 z ~ 8 objects identi- 
fied bv IBouwens et al.) (|2009af ). These galaxies are de- 
tected in the first epoch ultra-deep WFC3/IR imaging 
of the HUDF09 survey, which reaches ~ 29 mag (5it) in 
yi05>^i25-ffi6opfiotometry over an area of 4.7 arcmin^ cov- 
ering the Hubble ultra-deep field (HUDF: IBeckwith et al.l 
I2006D . The 16 z ~ 7 galaxy candidates are selected based 
on their [z^zo-Yioz) vs. {Y w^-Ji2^) c olors using the Ly- 
man Break Technique (e.g. ISteidel et al. 1996) and have 
an expected redshift distribution z ~ 6.5 — 7.5 with a me- 
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Fig. 1. — Surface brightness contours of the five brightest galax- 
ies in our sample and of a J125 stack of the remaining 11 fainter 
galaxies (last panel in lower right). The first five images are super- 
positions of lios, Ji25i and HiQo exposures, with the contour lines 
corresponding to /^yju = 23.5 — 25.5 mag/arcsec^ in steps of 0.5 
mag/arcsec^. The bar in the left corner indicates 2 kpc (physical) 
at z = 6.8 the expected mean redshift of these galaxies. AH images 
are I'.'S on a side. The size (FWHM) of the J125 PSF is shown as 
an inset in the lower right panel for comparison. 

dian at (z) = 6.8. Similarly, the z ^ 8 sample is selected 
based on their (Yios — J125) vs. (J125 — -ffieo) colors. 
For more infor mation on the survey and t he galaxy sam- 
ple we refer to lOesch et al.l ()2009b[ ) and iBouwens et all 
()2009ar ). We also include galaxies from 2; ~ 4 — 6 iden- 
tified as dropout galaxies in the optical HUDF data for 
comparison with the 2; > 6.5 population. 

We adopt = 0.3, = 0.7, Hq = TOkms-^Mpc-i, 
i.e. h = 0-7. Mag nitudes are given in the AB system 
(|Oke fc Gunnl[l983l) . We express galaxy UV luminosities 
in units of the charact eristic luminosity a t z ~ 3 being 
MiGooiz = 3) = -21.0 (iSteidel et al.lll999l) . 

2. MORPHOLOGIES AT Z - 7 

The 2 ~ 7 galaxies are extremely compact as can be 
seen from the contour plots in Figure [U where we show 
the individual summed Yios , J125 , and iJieo observations 
of the five brightest sources and additionally a J125 stack 
of the remaini ng 11 fainter ones ( for stamps of individ- 
ual sources see lOesch et al.|[2009bl) . As can be seen, the 
average 2; ~ 7 galaxy appears to be very symmetric and 
compact; with two exceptions, no extended features can 



be identified. The two exceptions are: 

a) The galaxy UDFz-42566566 is the brightest galaxy 
in our sample and consists of two clearly distinct com- 
ponents, separated by 2 kpc. These two each contribute 
about the same amount of light (1:1.2) with very similar 
colors. They have individual half-light radii of 0.5 kpc 
and 0.8 kpc, respectively, very similar to the compact 
galaxies in our sample. One interpretation for the origin 
of the individual components is therefore that they are 
in a merging phase. The linear geometry of the whole 
galaxy, however, may also suggest that the individual 
clumps are star-forming regions within a disk structure, 
similar to what has been found at z ~ 2 in observations 
and simulations. 

b) The galaxy UDFz-39557176 also consists of at least 
two components. The total light of this galaxy is dom- 
inated by a slightly elongated central structure, which 
has a fainter counterpart about 2 kpc away to the NW. 
The flux ratio of these two is 1:1.4. It is worth noting 
that this gala xy has been split into two sources in the 
iMcLure et all ([2009) catalog. However, the two knots 
are most probably physically connected and are about 
to merge with each other. The fainter component to the 
NW shows a significantly redder Y105 — J125 color com- 
pared to the central core by 0.2 mag, and also a redder 
J125 — HiQQ color by 0.1 mag. This may indicate that the 
second component consists of older stellar populations. 
However, a more speculative explanation could be that 
the second component is reddened due to dust from the 
central core. This would imply that these galaxies con- 
tain a more extended gas disk than what can be seen 
from their UV light. 

3. COMPARISON TO LEGS AT Z - 4 - 6 

In order to quantify the evolution of galaxy structures 
across cosmic time, we compare the 2; ~ 7 galaxies with 
LBGs identified at z ~ 4 — 6. We focus on three main 
aspects: (1) the size evolution, (2) the evolution of the 
average galaxy light profile, and (3) the evolution of the 
surface density of star-formation in these galaxies. 

3.1. Size evolution 

Galaxy sizes are measured using circular apertures 
containing 50% of the galaxies' light. We use the ob- 
served half-light radius from SExtractor, se' ^^'^ 
correct it for PSF broadening according to ri/2,sE = 
y^(''i/2 s_b)^ ~ ^PSF- "^^^ radius rpsF of a point source 

is 0'.'12 in the J125, and O'.'ll in the Y105 observations. 
These measurements are checked against the higher res- 
olution optical data for the 2: 5 population where size 
measurement from both ACS (1775) and WFC3/IR (Y105) 
are available. No significant differences are found, show- 
ing the robustness of these size estimates even for such 
faint, small galaxies and the validity of the simple PSF 
correction. 

As a seco nd check, we adop t size measurements based 
on galf it (|Peng et al.ll2002f l. Single sersic profiles are 
fitted to these galaxies with sersic indices fixed a.t n — 
1.5, the value measured for a mean stack of z ~ 4 LBGs. 
The half-light radii estimates of the best-fit models are in 
good agreement with the SExtractor measurements with 
a dispersion of (7^ ~ 0'.'05 and no bias. Similar results are 
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Fig. 2. — Top: Luminosity-Size relation for 2 ~ 5 — 8 galaxies 
in the WFC3/IR observations. The colored points are the z ~ 7 
(left) and 2 ~ 8 galaxies (right). Dark gray dots represent sam- 
ples of z ~ 5 — 6 based on WFC3/IR observations, while light 
gray squares are 2-^4 galaxies based on optical ACS data for 
comparison. All sizes are PSF corrected. The dashed line indi- 
cates the 50% completeness for galaxies with exponential profiles 
shifted to the absolute magnitude at the given redshift. The ver- 
tical dotted lines correspond to the edges of the luminosity bins 
with 0.12, 0.3, 1L*_3. The blue points mark the mean of the 2 ~ 7 
population in the given luminosity bin. Bottom: Comparison of 
the size distribution of 2 ~ 4 — 7 galaxies in the luminosity range 
(0.3-l)L*_3. The filled yellow histogram is for 2 ~ 7 galaxies, the 
others are for lower redshift LBGs (red: 2 ~ 6, green: 2 ~ 5, black: 
2 ~ 4). In order to increase statistics, this plot is based on GOODS 
for the galaxies of 2 ~ 4 — 6, where the sizes are measured from 
SExtractor in the 1775 (2 ~ 4) and 2350(2 ~ 5 — 6). The 2 ~ 7 
measurements are based on the J125. The main difference is the 
tail towards larger sizes at later cosmic times, while the peak of 
the size distribution changes little. 



found when using an average Sersic index n = 1, or n = 
3. However, galf it fails to return reliable measurements 
when the light profile is not well approximated by a single 
component fit and we use the SExtractor measurements 
as our fiducial ones. When appropriate, we also comment 
on the implications of using the galf it measurements. 

Due to selection effects, large, low-surface brightness 
galaxies will be missed in our catalog. We estimate this 
bias by inserting artificial galaxies of fixed profiles into 
the science images and rerunning the detection algorithm 
with the same parameters as for the creation of the origi- 
nal catalogs (for more details see lOesch et aLH2007l ). The 
galaxy profiles are chosen to follow an exponential Sersic 
function (n = 1) and a de Vaucouleur profile {n — 4). 

In the upper panel of Figure [2] we plot the half-light 
radii against the absolute magnitudes for the 2; 7 — 8 
galaxies and both z ~ 5 — 6 and z ~ 4 samples. As can be 
seen, the z 7 sample is complete for all galaxies at 1 kpc 
down to M1600 = —18.7, corresponding to Q.12L*^^. At 
these same luminosities the galaxies in the z ~ 8 sample 
are close to the completeness limit. We therefore do not 
include the z ^ 8 sample when fitting scaling relations 
as a function of redshift in our analysis. 

As can be seen from the figure, at fixed luminosities, 
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Fig. 3. — The evolution of the mean galaxy size across the 
redshift range 2 ~ 2 — 8 in two different luminosity ranges (0.3- 
1)L*=3. (top) and (0.12-0.3)L*^3 (bottom). The data from this 
work is plotted as red dots and is based on the PSF corrected SEx- 
tractor half-light radi i. Also reported are the measurements from 
IBouwens et al.l II2004I ) based on optical ACS data of the same lu- 
minosity range. These points are repeated in the lower panel for 
comparison. Mean half-light radius measurements from Sersic fits 
with galfit arc shown as orange circles, offset from their redshift 
by 0.1. The dashed lines indicate the scaling expected for a fixed 
halo mass (oc H{z)~'^^^ ; black) and at fixed halo circular velocity 
(oc H{z)~^; gray). The red lines correspond to the best-fit to the 
observed evolution described by (1 -|- 2)"*", with m = 1.12 it 0.17 
for the brighter luminosity bin, and m = 1.32 it 0.52 at fainter 
luminosities. These are formally identical and are consistent with 
the m ^ 1 value derived previously. The short black line at 2 ~ 7 
indicates the mean sizes found when projecting the observed 2 ~ 4 
sample to 2 ~ 7, while keeping their luminosities and physical sizes 
constant in order to test for measurement biases due to cosmolog- 
ical surface brightness dimming. 

the sizes evolve only little from z ~ 4 — 8. However, 
the lower redshift population contains an extended tail 
towards galaxies with sizes > 1.2kpc, which is not seen 
at higher redshifts as shown in the lower panel of Figure 
[21 In order to increase the statistics, the z ~ 4 — 6 
galaxy sampl es are taken frorn GOO DS based on ACS 
imaging (e.g. iGiavalisco et al.l l2004allbl V The peaks of 
the distributions are between 0.8 — Ikpc for all redshifts. 
The z 4 population contains a significant population of 
galaxies with > 1.2kpc, which are completely absent at 
z ~ 7. Firm confirmation of the lack of such large objects 
at z 7 will require larger area surveys, but the trend 
seen from z ~ 6 to lower redshift is very suggestive. The 
above indicates a scenario in which these primordial star- 
forming galaxies are formed as small clumps with a size of 
^ 0.4 — O.Skpc independent of luminosity. Towards later 
cosmic times hierarchical build-up or enhanced accretion 
of cold gas starts to produce larger galaxies. 

The evolution of the average sizes of galaxies from z ~ 
2 — 8 are plotted in Figure [3] for two different luminosity 
ranges, (0.3-l)L*^3 (top) and (0.12-0.3) L*^3 (bottom). 
The expected size scaling from semi-analytical models 
is H{z)~^ , at fixed halo circular velocity, or H{z)~'^/'^ at 
fixed halo mass (Mo et al. 1998), where H{z) is the Hub- 
ble parameter at redshift z which scales as ~ (1 -t z)'^/^ 
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Fig. 4. — Apparent intensity profiles of the stacked z ^ A — 7 
galaxies in two different luminosity ranges, (0.3-l)L*_3 (left) and 
(0.12-0.3)Z;*_3 (right). The blue points are the data, while the 
black line is the PSF from stars in the J125 observation. The best- 
fit model of the z ~ 7 population is shown as a solid green line. 
Models for lower redshift galaxies are shown as dashed lines in 
red (z ~ 6), cyan (z ~ 5) and black (z ~ 4). For galaxies with 
luminosities (0.3-l)L*_3, the z ~ 4 — 5 population shows clearly 
extended wings with respect to higher redshift galaxies, while in 
the lower luminosity bin, the evolution is slower. Note that this 
is consistent with increasing physical sizes towards lower redshifts 
due to increasing angular diameter distances. 

at z > 2. The observed size evolution at fixed lumi- 
nosity is fitted with a scaling of the form (1 + z)~™. 
The two fits are formally identical with m = 1.12 ± 0.17 
and m — 1.32 ± 0.52, respectively. This is in agree- 
ment with previous estimates where the s izes were found 
to sc a le rou ghly according to (1 + z)~^ (iBouwens et alj 
2OOI [2OOI How ever, H{z)-^ (jFerguson et al.l 120041 : 
Hathi et "aLll2008bl ) cannot be ruled out as both scalings 
are very flat over our redshift range and diverge only at 
lower z. 

Note that the observed evolution in sizes is not an ar- 
tifact of cosmological surface brightness dimming. We 
verified this by artificially redshifting our z ~ 4 sample 
(1775 band) to z ~ 7 ( J125 band) and remeasuring their 
sizes with SExtractor. For the more luminous galaxies, 
the input size is recovered perfectly, while for the lower 
luminosity sample only a small bias of ~ 10% towards 
smaller sizes is found due to the flux loss in the galaxy 
wings. 

3.2. Average Profiles 

In order to test the uniformity of average galaxy sizes 
over these redshifts, we measure their average profiles 
based on stacking galaxies from fixed luminosity ranges 
and fitting these profiles with galfit. For the z ^ 7 
galaxies this is done in the J125 band, corresponding 
closely to restframe 1600 A. The galaxy stamps are 
cleaned o f close neighbors us ing the cleaning algorithm 
of ZEST+ (ICaroUo et al.ll2009L in preparation). The algo- 
rithm masks neighboring sources in the images and re- 
places contaminated pixels with background noise. The 
point-spread function (PSF) was estimated on three stars 
which showed no saturation in individual exposures. We 
have also ensured that the central pixels of the stars have 
not been affected by cosmic-ray rejection in the data re- 
duction pipeline. 

The above procedure is repeated in the Y105 band for 
samples of z ~ 5 and z ~ 6 galaxies, extracted from the 
optical HUDF catalogs. The Y105 corresponds to 1750 A 
and I5OOA restframe at z ~ 5 and z ~ 6 respectively. 
For the z ~ 4 population the 1775 band corresponds most 



closely to the I6OOA wavelength. Therefore morphologi- 
cal k-correction effects are minimal. 

A Sersic profile is fitted to these stacked images from all 
three redshift bins and both luminosity ranges. Through 
simulations, we have verified that the stacking procedure 
reproduces the av erage profiles of a ga laxy population ac- 
curately (see also lHathi et al.|[2008al) . In order to com- 
pare the different profiles directly, we have convolved the 
best-fit models with the J125 band PSF and show the 
result in Figure 2) The faint z ~ 7 galaxies are only 
marginally resolved. The average profile is remarkably 
similar at all redshifts from z ^ 4 to z 7, except for 
the appearance of extended wings at z ~ 4 for all lu- 
minosities, and extended wings at z ~ 5 for the more 
luminous galaxies. But even the similar angular size in- 
dicates increasing physical sizes towards lower redshifts 
due to increasing angular diameter distances (30% larger 
at z ~ 4 vs. z ^ 7). 

The steepness of all these profiles suggests that star- 
formation is largely centrally concentrated in these galax- 
ies and not widely distributed in a large, clumpy disk. 

3.3. SFR surface density 

Given the slow evolution of the sizes of LBGs, it is 
interesting to estimate the change in the average sur- 
face density of star-formation, "Esfr- The upper panel 
of Figure [5] shows this as a function of absolute mag- 
nitude. The star-formation rate (SFR) has been esti- 
mated based o n a simple convers ion of the UV lumi- 
nosity to SFR (Ma dau et a"l]|1998f) assuming a Salpeter 
initial-mass function. For the z ~ 7 galaxies the SFRs 
range from ^ 1 to lOM0yr~^. Furthermore, our results 
suggest that the average "Ssfr remains relatively con- 
stant for the whole redshift range from z '--^ 7 to z ~ 4. 
This is shown for galaxies with luminosities in the range 
(0.3-l)i*^3 in the lower panel of Figure El Filled red 
circles represent the actual measured values, while the 
orange open circles are corrected for du st extinction us- 
ing the formula bv lMeurer et al.l (jl999l ) and UV contin- 
uum slopes measured in this same luminosity and red- 
shift ranges from Bouwcns et al. (2009b). In agreement 
with the size evolution, there is a small trend towards 
lower apparent ^sfr at lower redshifts. However, when 
corrected for dust extinction, the star-formation surface 
density is consistent with being constant over the entire 
redshift range z ^ 3 — 7. 

A possible explanation for the constancy of the star- 
formation surface density is that the average star- 
formation efficiency is very similar in all these galax- 
ies and that feedback effects change the mode of star- 
formation only little in these primordial galaxies. Note 
that this result is consistent with a constant limiting 
surf ace brightness frorn z ~ 7 down to lo cal galaxies 
(see iMeurer et all Il997t iHathi et all l2008b[ ) . However, 
there are exceptions from this typical mode of star- 
formation. For example, the average surface density 
of star-formation is almost three orders of magnitudes 
larger in hyper- starbursts in quasar hosts at similar red- 
shifts (z = 6.42: iWalter eTaH l2009l ) . 

4. SUMMARY AND CONCLUSIONS 

We have used the ultra-deep WFC3/IR observations 
of the HUDF09 program to study the structures and 
morphologies of z > 6.5 galaxies previously presented in 
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Fig. 5. — Top: The surface density of star-formation against 
luminosity. As in Figure [2] the gray background points represent 
the 2: ~ 4 — 6 population in the WFC3/IR data and the colored 
points are the 2 ~ 7 (left) and 8 population (right). Bottom: 
The mean surface density of star-formation rate against redshift. 
Red dots correspond to the actual measurements, while the orange 
circles are corrected for dust absorption (Bouwcns ct al. 2009b). 
The black line is the average for the uncorrected values and the 
dashed line the same for the dust corrected ones. In both cases, 
the average Tigpn is consistent with being constant over the entire 
redshift range 2 ~ 4 — 7. 



lOesch et all (|2009bf ) and lBouwens et all (|2009a| ). These 
galaxies are all extremely compact, with an average size 
of 0.7 ± 0.3 kpc. Only two out of the 16 z ^ 7 galax- 
ies show a double core. This fraction is slightly smaller 
than the ~ 30% of LBGs with disturbed m orphologies 
and d ouble cores found at lower rcds hifts (e.g. iLotz et all 
120061: IConselice fc Arnoldl I2009t iPettv et all I2009D . al- 
though they are still consistent within the small number 



statistics. The z ~ 7 galaxies are extremely similar to the 
z ~ 6 population, both in sizes and in their average light 
profiles, showing that galaxy evolution proceeds slowly 
over the ~ 170 Myr from z ^ 7 — 6. 

By comparison to LBGs down to z ^ 4 only a very 
slow size evolution is found, following (1 -I- z)"™, with 
m = 1.12 ± 0.17 for galaxies of luminosities (0.3-l)L*^3. 
Fainter galaxies down to 0.12L*^3 follow a similar scal- 
ing with m = 1.32 ± 0.52. Additionally, the mean star- 
formation surface densities of LBGs are found to be con- 
stant over the entire redshift range z ~ 4 — 7, which may 
be explained by largely constant star-formation efficien- 
cies at these early epochs. 

Similar exponents for the size scalings with redshift 
have been fouiid for d isk galaxies between z ~ — 3 (e.g. 
iBuitrago et al.l l2008l ). as well as from semi- analytical 
modeling including concentrated dark-matter halo pro- 
files (e.g. ISomerville et al.|[200l iFirmani et al.ll2009[ ) 

By extending the present study to larger samples of 
z ~ 7 galaxies which will become available in the near 
future it will be possible to put more stringent constraints 
on galaxy evolution models and to shed more light on the 
question of the main driver of star-formation in LBGs at 
these early epochs. 
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